Influence of exciton-phonon interaction on long energy transport in J - aggregates. 
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This paper presents a theoretical model intended to address the question of energy transfer in 
two-dimensional molecular assemblies such as Scheibe aggregates. A new phonon-exciton interaction 
is introduced to explain the exciton width in J aggregates. It is shown that the long range energy 
transport can occur for weakly interacting acceptors. 



I. INTRODUCTION 

Compact and regular two dimensional arrangements of oxycyanine dye bound to hydrophobic carbon tails show 
remarkably narrow absorption and fluorescence peaks [1-5]. In these so called J-aggregates it is believed that the 
excitation of a dye molecule becomes delocalized (exciton) due to the strong dipole-dipole interactions between 
delocalized electrons in adjacent dye molecules. It was also shown [3-5] that these structures can support practically 
loss-free energy transport over large distances. Molecules of a suitable acceptor (a molecule similar to that used 
in J-aggregate, but with thiocyanine instead of oxycyanine) were able to absorb almost all the energy, even at 
concentrations as low as 1:10 000. It is, however, astonishing, that energy transport is more efficient when the 
coupling to the acceptor is weak [5]. Quantum mechanical mechanism of this effect was developed in [7,8] in terms of 
exciton trapping by shallow bound state of large spatial extend. However, the theory assumed 'ideal' excitons, which 
would imply extremely small absorption width in J-aggregate. Clearly a dissipation mechanism is needed to explain 
the observed experimentally absorption width, which, moreover, should not destroy the delicate exciton trapping 
needed to explain long range energy transport. Investigation of possible dissipation mechanisms is the aim of this 
paper. It will be shown that the usually assumed [6] exciton-phonon coupling is not up to the task and that only 
a mechanism destroying excitons will do. We propose, therefore, a new phonon-exciton interaction, which has the 
required properties and, moreover, leads to a solvable model. Then we demonstrate that our interaction permits a 
long range energy transfer in two dimensional structures. 

II. IDEAL EXCITONS COUPLED TO PHOTONS 

We start with the one dimensional Hamiltonian for the excitons [6,7]. We have 

oo 

H= J2 pA+A n -JA+(A n . 1 + A n+1 )}, (1) 

n— — oo 

where A^, A n are the creation and annihilation operators, respectively, of an excitation at site n, is the excitation 
energy of a single dye molecule and J is the resonant interaction energy between adjacent sites. We use one dimensional 
formalism for simplicity's sake, two dimensional generalization is immediate. Hamiltonian (1) is diagonalized in 
momentum space. Introducing Bloch waves with quasi momentum k we have 

A n = -1= / dke ikn A k , (2) 



/2tt 

A »= E e "*" A »- ( 3 ) 



Our Hamiltonian (1) is now 
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H= f n(k)A+A k dk (4) 

J — 71 

with the dispersion relation 

Q(k) = Q - 2Jcos(fc). (5) 

Again for simplicity we use scalar photons with the coupling to the dye within rotating wave approximation [8], so 
the Hamiltonian for the exciton-photon system is 

H = J dkfl(k)A+A k + J d 3 pe( P )c+c p + (6) 
d 3 pdkg(p)S( Pll - k){A k c+ + c p A+), 



where the second integral is the energy of the photons with the usual dispersion relation 

e( P ) = c\p\ (7) 

and the last term is the photon-exciton interaction. Note that the sum of exciton and phonon numbers is conserved (by 
the virtue of rotating wave approximation) and that we have delta of momentum conservation along the chain. At the 
momenta of interest the molecular formfactor g(p) is practically a constant, which we denote go, and its momentum 
dependence is of importance in loop integrals only: it gives the natural line width of an atomic state, much smaller 
than the experimentally observed width [8]. It is important to realize that since the wavelength of visible light is of 
the order 1000 A, this corresponds to momenta of the order 10~ 3 A -1 , so to a good approximation light produces 
excitons with negligible momentum, i.e. from the bottom of the band. In this form the model is analytically solvable. 
Working within one photon-one exciton sector we use for the wave function the ansatz 

|V >= f dk<p(k)A+\Q >+J d 3 P ^ Pa (p)c+\0 >, (8) 

where |0 > is the Fock vacuum (which is an exact eigenstate of the theory). This leads to a system of equations for 
the unknown functions ip(k) and tp Po (p) 

zip(k) = Q(k)ip(k) + J d 3 pg(p)S(p\\ - k)ip Po (p) , Q , 

z ^po(p) = £(p)ip P0 (p)+g(pMp\\)- 

where the photon wavefunction ip Po (p) is a sum of plane and scattered wave 

Tp Po (p) = 5 3 (p- Po ) + *+(p), (10) 

where po is the momentum of the incoming light. System (9) is easily solved and we get for the experimentally 
observed photon wavefunction 

s (pq\\-p\\) g(p)g(po) 

(e(po)-e(p) (z-n( Pll )-G( Pll )y 

As we see, the scattered wave has a cylindrical outgoing wave (note the one dimensional delta function) multiplied 
by scattering amplitude. The absorption intensity is proportional therefore to 

g{p)g(po) 



(p) = ppLJ^L (ii) 



z - Sl(j>\\) - G(p\\) 



(12) 



where we have introduced the notation 

„2 



G(k) = [ d ^^~ k \ (13) 
J z- e{p) - is 

for small k (our case) the real part of G(k) is negligible, while the imaginary part is 

G(0) = tirgl (14) 
The absorption probability has therefore the familiar lorentzian shape 

1 



(z-m) 2 + K 2 g 4 ' 



(15) 



Due to the smallness of go [8], the absorption peak has a finite, albeit very small with. In order to get a more 
realistic value, we need a dissipation mechanism, i.e. interaction with phonons. 
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III. PHONON-EXCITON INTERACTION 



Like everybody else we assume that excitons interact with phonons, since the J-aggregate structure is not very stiff. 
We take the phonon Hamiltonian in the familiar form 



H = / uj(q)a q a q dq. (16) 

J —IT 

where a+, a q are the (bosonic) creation and annihilation operators for phonons with the dispersion relation 



u(q) = v /2w (l - cos(g)) ^ s/uo\q\, (17) 
One usually assumes a trilinear coupling exciton-phonon [6] 

Hi = J X (k, q)A+ +q A k (a q + at q )dkdq. (18) 

However, this coupling does not change the exciton number: This interaction can only push excitons to the bottom 
of the band, but never destroy them! Since excitons are created with negligible momentum anyway, this mechanism, 
although certainly present, cannot contribute significantly to the observed light absorption width. We need therefore 
an interaction which will explicitly destroy excitons, regardless their energy. The simplest possibility is the annihilation 
of an exciton into 2 phonons. After all, the excited dye molecule contains a dclocalizcd electron, which can not only 
reradiate light, but also give up its energy to the vibrations of the whole structure. We of course realize that it 
is difficult for two phonons to carry away 2-3 eV of energy, so we call them 'effective'. We therefore propose the 
interaction in the form 

Ht = J x(fc, q)(A k a^_ q a+ + A\ a k - q a q )dkdq. (19) 

where x(^? <?) IS t ne formfactor, which, for the sake of simplicity, we assume is constant. Our interaction has the great 
advantage that it leads to a solvable model. Our excitons coupled to photons and phonons are described by a free 
part 

H = J tt{k)A+A k dk + j io{q)a+a q dq + j e(p)c+c p d 3 p (20) 



and the interactions 



H i = J x(A k a^_ q aj + A+a k - q a q )dkdq + 

J g P {A k c+ +c p A+)6{k-p ll )dkd 3 P . (21) 
One exciton sector is solved by the ansatz 

|*> = J dk(p(k)A+\0> + J dq 1 dq 2 ^p h (q 1 ,q2)a+a+\0> + 

J d 3 p* po (p)4\0> (22) 
which leads to a coupled set of equations 



f zip(k) = n(k)ip(k) +2J dq X (q)ip P h(q, k - q) + 

f d 3 pg p ii> P0 (p)S(pu - k) 

ztppoip) = £(p)V>po(p) +9 P <P(P\\) 
I zip ph (q 1 ,q 2 ) = X(q2)<p(qi+q2) + {u(qi)+u(q2))ip p h(qi,q2)- 



(23) 



where the energy is 

z - e(po) (24) 



3 



The scattered photon wave function can be written in the form 



y+ ( v ) = g Po g p S(Po\\-P\\) (25] 

?° [P> { z -Q(p ll )-r z (p ll )-G(p ll )}{e(p Q )-e(p)y { > 



where we have introduced the function 

X 2 (q) 



/7T 
■ 
-7T 



Z — Lo(q) — Lo(k — q) — ie 
Absorption cross section is therefore proportional to 

1 



dq. (26) 



z-n(k)-T x (k) ' 

Remembering that excitons are created with zero momentum we need the function (26) at k=0, we have 



(27) 



r,(o) 



= 2x 2 f ■ ■ (28) 

J-t z-2u{q) - ie 



Assuming linear dispersion for the phonons (cf. eq. 16) we can perform the integration explicitly 

1 



(--^(0) + ^ln|l-^o|) 2 + ^' 



(29) 



The logarithm does not play, numerically speaking, any significant role, therefore we have obtained a familiar 
lorentzian shape. Its width is given by 

/m 2 r,(0) = (30) 

and depends on \ an d w o- Since experimentally it is of the order of 0.026 eV, we can estimate \ if we assume iv 0} 
which we arbitrally have taken to be 4.0 eV. The O(0) parameter gives the position of the absorption peak, taking it 
to be 3.104 eV we get an excellent agreement with the experiment. 

IV. ACCEPTOR IN THE PRESENCE OF PHONONS 

Presence of an acceptor molecule, which we place at site n = 0, introduces an additional term in the exciton 
Hamiltonian, which now reads 

oo 

H= m+A n -JA+(A n _ 1 + A n+1 )]-(n-Cl Q )A+A , (31) 

n=— oo 

where f2 is the acceptor excitation energy [7] . In momentum space this is equivalent to the introduction of a separable 
potential with constant formfactor f(k) 

H = J n(k)A+A k dk + X J dkf(k)A+ J dkf{k)A~ k (32) 

and 

A = (33, 

As shown in [7,8] for small values of lambda, leading to a shallow bound state, we get large absorption, while 
paradoxically large lambda does induce strong absorption. The question now is: will the interaction with phonons 
destroy this delicate mechanism? A satisfying feature of our interaction is that the model is still analytically solvable. 
The full interaction is now 
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Hi = I x{q){Akat_ q a+ + A^a k ^ q a q )dkdq + 
g p (A k c+ + c p A+)5(k - Pl] )dkd 3 p + 
A fdkf(k)A+ fd~kf(~k)A~ k , 



(34) 



which gives now a set of equations in one exciton - one photon - two phonons sectors (which are coupled) 

' zip(k) = tt(k)<p{h) + 2j dqx(q)ip p h(q, k - q) + 

\f(k)jd~kfCkM~k)+Jd 3 P gM P )S( Pll -k) 

ZTppo(p) = e(p)i>p (p) + 9p<P(p\\) 
^ zip ph (qi,q 2 ) = x{ii)^p{qi +92) + (Mgi) + u{q2))i) P h{<ii,q2)- 

Again the system can be solved and we get the scattering part of the photon wave function 

9 Po 9p5(Po\\ ~P\\) 



(35) 



{z - n(p,|) - r,(p,|) - G?(p||)}{e(po) - e(p)} 
9p9 P0 f 2 {P\\) 



(36) 



{ e (po)-e(p)}{^-^0||)-r ;e (p||)-G(p||)}{^-0(po ll )-r z (po lt )-G r Oo ll )} 

dkp(k) 



I 



n(k)-r z (k)-G(k) 



(37) 



showing the effects of acceptor and phonon interactions in all its gory details. The last line is the effect of the acceptor. 
We will obtain a significant enhancement if this term is small. As assumed above we can take 



T z (k) ~ r z (0) = iSexp 

where S exp is the experimental width. The whole behavior is therefore controlled by the integral 



so*) = r 

J —7T 



dk 



Z O(fc) i&exp 



(38) 



(39) 



For suitable values of coupling lambda real part of E can be canceled, leading to a large enhancement, provided 
that the imaginary part is reasonably small. Inspection of FIG. 1 
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FIG. 1. Function E(z), 1 dim case. 
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, where the y axis is at the threshold for exciton production, shows indeed that large values of j (i.e. lambda 
small!) are favored. Absorption spectra with j small (FIG. 2) and j (FIG. 3) show an enhancement of the order of 
50%. 
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FIG. 2. Absorption cross section - weakly bound acceptor (1 dim) 
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FIG. 3. Absorption cross section - strongly bound acceptor (1 dim) 



The effect is there, but disappointingly small. The reason lies in the large value of the imaginary part of sigma (cf. 
FIG. 1). We must notice however that this is a one dimensional model. In two dimensions we would get formally the 
same expressions, but all integrations are 2 dimensional. In particular the function sigma, after angular integration 
would be 

= f J, d , k c (40) 

J Z - Si(fc) - Iddos 

and has integration measure kdk, which suppresses small k - just what is needed. The imaginary part at threshold is 
now strongly reduced (FIG. 4). 
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FIG. 4. Function £(z), 2 dim case. 



The absorption enhancement at small lambda now is of the order of 700% (FIG. 5,6). 
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FIG. 5. Absorption cross section - weakly bound acceptor (2 dim) 
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FIG. 6. Absorption cross section - strongly bound acceptor (2 dim) 
It is seen that the fact that J aggregate is a two dimensional structure is of paramount importance. 
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V. CONCLUSIONS AND OUTLOOK 



We have shown that our postulated mechanism: exciton annihilation into 2 phonons leads to a solvable model. 
It describes correctly the width of J-aggregate absorption spectra. In the case of introduced acceptors, the large 
absorption on shallow bound states is confirmed, but two dimensional geometry seems essential. Of course full two 
dimensional calculation is needed and will be presented shortly. In order to describe fluorescence spectra, we have 
to include the usual energy loss by phonon radiation mechanism [6]. In conclusion let us remark that absorption 
energy on shallow bound states, postulated already in [7] not only survives the introduction of realistic dissipation 
mechanisms, but acquires new features allowing for a better agreement with experiment. 
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